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=^ : ^he re^„t reports of temporal and spectr-al pecuHanttes in the early stages of some after- 

glows suggest that we may be wrong in postulating a central engine which becomes dormant 
' after the burst itself. A continually decreasing postburst relativistic outflow, such as put out 

by a decaying magnetar, may continue to be emitted for periods of days or longer, and we 
(N ■ ai-gue that it can be efficiendy reprocessed by the ambient soft photon field radiation. Photons 

produced either by the postexplosion expansion of the progenitor stellar envelope or by a bi- 
nary companion provide ample targets for the relativistic outflow to interact and produce high 
, energy 7-rays. The resultant signal may yield luminosities high enough to be detected with 

' the recently launched Integral and the Glast experiment now under construction. Its detec- 

, tion will surely offer important clues for identifying the nature of the progenitor and possibly 

constraining whether some route other than single star evolution is involved in producing a 
rapidly rotating helium core which in turn, at collapse, triggers a burst. 
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1 INTRODUCTION 

The typical GRB model assumes that the energy input episode is brief, typically tgrb < 1 — lO'^s (see Meszaros 2002 for a recent review). 
However, peculiarities in the early stages of some afterglows, e.g. GRB 021004, have served as motivation for considering a more extended 
input period in which the energy injection continues well beyond tgrb (Fox et al. 2003). The enduring activity could in principle emanate 
from the sluggish drain of orbiting matter into a newly formed black hole (e.g. MacFadyen & Woosley 1999) or from a spin down millisecond 
super-pulsar (e.g. Usov 1994), which could produce a luminosity that was, still one day after the burst, as high as L '--^ 10*^ erg s~^. Collapsar 
(Woosley 1993, Paczynski 1998; MacFadyen & Woosley 1999; Aloy et al. 2000) or magnetar-like GRB models (Usov 1994; Thompson 1994, 
Wheeler et al. 2000) provide a natural scenario for a sudden burst succeeded by a more slowly decaying energy release (Rees & Meszaros 
2000; Ramirez-Ruiz, Celotti & Rees 2002). 

The power output would be primarily in a magnetically driven relativistic wind, which would be hugely super-Eddington during the 
time scales discussed here. Its luminosity may not dominate the afterglow continuum (Dai & Lu 1998; Zhang & Meszaros 2001), but we 
argue that it could be efficiently reprocessed into 7-rays. This could be due to the interaction of the postburst relativistic outflow with the 
dense soft photon bath arising either from a stellar companion or from the postexplosion expansion of the remnant shell or supernova (Hjorth 
et al. 2003; Stanek et al. 2003). The detection of such scattered hard 7-ray radiation would offer the possibility of diagnosing the nature of 
the precursor star and the compact object that triggers the burst. 



2 BINARY SYSTEM WITH A DECAYING MAGNETAR 

In the generic pulsar model the field is assumed to maintain a steady value, and the luminosity declines as the spin rate slows down. However, 
during the early stages, the magnetic field strength might decline more rapidly than the slowing down timescale. The power output in this 
case declines in proportion to B^. The spindown law is given by —Ifltl = _B^i?®r2*/(6c'^) + 32G/^e^f2®/(5c^) (Shapiro & Teukolsky 
1983), where Q is the angular frequency, B is the dipolar field strength at the poles, R is the radius of the light cylinder, / is moment of 
inertia, and e is ellipticity of the neutron star. The above decay solution includes both electromagnetic (EM) and gravitational wave (GW) 
losses. At various times the spin-down will be dominated by one of the loss terms, and one can get approximate solutions. When EM dipolar 
radiation losses dominate the spin-down, we have = f2o(l + t/tm)~^^'^ ■ Here 
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Figure 1. The luminosity emitted by the (postburst) relativistic outflowing wind through the Compton-drag process. The power output of the magnetically 
driven relativistic wind declines as the spin rate slows down (see eq.|2) and B is assumed to maintain a steady value (dashed line). Lni,o — 10*^ erg s~^, 
tni — 10^ s, and ^ = 1. The luminosity of the 7-ray emission generated as the relativistic outflow interacts with the soft photon field of the companion is 
shown for ?=1,10 and for two different values of the inclination angle (i.e. i). The period of the binary in the circular orbit is Tp = 12 hours. The emitted 
luminosity was derived under the assumption that the relativistic wind is both steady and mono-energetic. In reality, the properties of the continuing power 
output could be more complicated. Such effects will certainly enrich the temporal and spectral dependence of the resultant 7-ray emission. 
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is the characteristic time scale for dipolar spindown (e.g Dai & Lu 1998), B15 = B/(10^^G), Pq is the initial rotation period in milliseconds. 
When GW radiation losses dominate the spin-down, the evolution is f2 ~ f2o(l + ^/igw)^^''*, where tg„ ~ Is/^^^Pq _3(e/0.1)~^. GW 



spindown is important only when the neutron star is born with an initial fio > ^* 
the continuous injection luminosity is given by 

I/m(t) = Lm,0(l +^/^m)"^ 



lO^'s"^ (e.g. Blackman & Yi 1998). When < n.. 
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If the neutron star is bom with flo > fi*, the timescale for GW-dominated regime is short so that Q. will be damped to below Jl, promptly 
in a time tt = [(f2o/f2*)* — l]igw After Q, < Q.,, GW losses decrease sharply, and the spin-down becomes dominated by the EM losses. 
The injection luminosity can therefore be divided into two phases, i.e. L — I/m,o/(l + i/igw) for t < t* or L = Lni,*/[1 + (i — t*)/tm,*^ 
otherwise, where Lm,, = 10.1/ {2t^) and tm,* = 2,c^ I / {B"^ R^Q.,). 

The bulk of the magnetar energy Lm would be primarily in the form of a magnetically driven, highly-relativistic wind consisting of e^, 
e'^ and probably heavy ions with Ln ~ (^Lm and C < 1- We envisage that the burst is triggered by the collapse of a massive star whose 
helium core is presumed to be kept rapidly rotating by spin-orbit tidal interactions with its binary companion. The relativistic (postburst) wind 
would then escape the compact remnant while interacting with the soft photon field of the companion with typical energy 9 = kT / {me(?) ~ 
10"'' — 10~^ (i.e. optical and UV frequencies). The scattered photons whose energy is boosted by the square of the bulk Lorentz factor of 
the magnetized wind (i.e. Q = 2^^d) propagate in a narrow 7"^ beam owing to relativistic aberration. In the case of a mono-energetic 
isotropic magnetar wind (which is likely to be undisturbed even in the presence of a strong mass outflow from the binary companion; i.e. 
L{i /c S> vM), and assuming that the relativistic particles and soft photons are emitted radially from their source, the energy extracted from a 
relativistic particle by the inverse Compton process in the Thompson regime is given by d7(r, ip) ~ — L,(tt7^ (1 — cos uj^'dr / {A-mns5'^(?'), 
where L, is the luminosity of the stellar companion, uo is the angle between the photon and particle directions before scattering, and r denotes 
the distance to the magnetar from the volume region where the 7-ray radiation is generated. Here 5 — \J A? sin^ ip + [r — A cos ipY defines 
the distance from this volume to the companion, A the binary separation, ip the angle between the lines connecting the binary members and 
the observer, and 5 cos a; = r — A cos p. The decrease of the particle's Lorentz factor from the magnetar to infinity is 



7 

where 



30 



A 

Rq 



(4) 



(5) 



Identifying young GRB fossils 3 



10" - 



oOMo 



10-6 



103 



w 10^ w 10^ 10^ 



Figure 2. Efficiency in extracting energy from the relativistic outflow as a function of tlie energy of the scattered 7-ray radiation. The properties of the observed 
signal are linked to both the evolutionary history of the companion and the bulk Lorentz factor of the relativistic wind (i.e. 7). The main sequence radii and 
luminosities of the companion star are calculated using the analytic function derived by Tout et al. (1996). We assume A = 3i?, . The dashed line denotes 
the threshold energy at which the beamed photons are absorbed by the soft field radiation. The two vertical solid lines mark the lower limits to the maximum 
redshift to which TeV 7-ray sources are visible, based on limits to the background density of IR/UV photons (Biller et al. 1998). This has been done for the 
7-ray energy range of 0.4-10 TeV, assuming that a source remains visible out to an optical depth of 2 and that the majority of the IR was produced prior to the 
epoch of the sources under study. 



denotes the efficiency in extracting energy from the relativistic outflow (mainly composed of e^), and 'I'(<^) = {3[tt — (/?] — i sin2</3 — 
4 sin (p)/{2 sin (p) is the beaming function of the generated 7-ray radiation^. The total luminosity emitted by the relativistic outflowing wind 
through the Compton-drag process in the direction of the observer is highly anisotropic: I/r(<^) = [A7(<^)/7]Ln. Lr strongly varies with 
ip which in turn changes periodically during orbital motion (i.e. cos p — sin i cos u for a circular orbit, where i is the inclination angle and 
u is the true anomaly). Fig. Q shows the luminosity carried by the scattered photons for various positions of the companion and different 
assumptions regarding the magnetar luminosity. 

The Compton drag process can be very efficient in extracting energy from the outflowing relativistic wind provided that ? > 1 (see 
equation|5j- ? is of course uncertain, but this number does not seem unreasonably high for a progenitor star with a massive binary companion, 
and suggests that our fiducial value of in for the overall luminosity need not be an overestimate (Fig.|2}. PSR B1259 - 63 is an example 
of a system where the mildly relativistic outflow from the aged pulsar is thought to interact with the soft photon field radiation produced 
by its high mass binary companion (Chernyakova & Illarionov 1999). Since ? could be less than 10""^ for a stellar companion with modest 
luminosity (or a binary with A R^,), this suggests that we cannot rule out the possibility that part of the 7-ray continuum could still come 
from the intrinsic pulsed emission from the magnetar itself. For a typical aged pulsar this mechanism yields Lr /Lq ~ 10~^ — 10"'^ (Arons 
1996). The time dependence and spectral properties, in this latter case, will be very distinctive and such effects should certainly be looked 
for. 



2.1 Spectral attributes 

If the scenario proposed here is indeed relevant to an understanding of the nature of GRB progenitors, then its existence becomes inextri- 
cably linked to both the evolutionary history of the companion and the characteristics of the hydromagnetic wind. The postburst outflow 
from the compact remnant (which we assume to be highly relativistic) propagates in the soft photon field of the stellar companion where 
U{r,zu) = mec'^t^ni^{r,zu) is the photon energy density at the location r. As the stellar companion emits a black body spectrum, of 
effective temperature 6, , the local photon energy density is given by 

where vj is the soft photon energy in units of m^f? . The scattered photons are boosted by the square of the Lorentz factor so that the local 
spectrum has a black body shape enhanced by 7^ . 

The luminosity of the scattered emission moving along a radial trajectory at an angle is Lr(7, V'; = ai^n^(\—fi co^Lo^a^K^dV , 
where e is the energy of the scattered photons in units of ra^f? , and ctkn is the Klein-Nishina cross-section (Jauch & Rohrlich 1976). As 

^ Owing to the fact that the optical star is not a point source, ^(}f>) is clearly only accurately for Rt/5 < ip < it — Rt/5, where ij, is the radius of the 
stellar companion. 
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Figure 3. Representative spectra produced by the interaction of tlie relativistic outflowing wind witli the radiation of the stellar companion. The power output 
of the magnetically driven relativistic wind declines as (1 + t/tm)"^, with Lm,o — 10^^ erg s" i, tm — 10'^ s, f = 1, and t = Tp = 12 hours (circular 
orbit). We assume z = 0.5 (Hq = 65 kms"i Mpc^^, = 0.3, and Oa = 0.7). The luminosity of the 7-ray emission is shown for 7 = 10* (dotted 
line), 10^ (solid line), 6», = 2.5 X 10"'', and f = 1. The thick grey lines are the predicted sensitivity of a number of operational and proposed instruments in 
7-ray astrophysics (Morselli 2002). 

can be seen in Fig.|3| the resulting spectrum is the convolution of all the locally emitted spectra (i.e. dn^ [r, -go]) and it is not one of a 
blackbody. Note that Klein-Nishina effects are important for incoming photon energies such that 1177(1 + cos ip) > 1. The maximum energy 
of the scattered photons in this regime is jmeC^. 

A further effect which may strongly affect the observed spectrum is the production of pairs through photon-photon collisions, 
pairs can be produced by scattered photons interacting with the isotropic companion emission or with each other. Photon collisions within the 
beam itself occur between photons of equal age and can only affect the high energy tail of the spectrum provided that 7S* > 1/3 (Svensson 
1987). Let us thus consider in turn the role of scattered and companion radiation as seed photons for this process. The interaction between the 
7-rays produced by the Compton-drag process and photons emitted by the companion star would occur at large angles, resulting in an average 
energy threshold of e^-, > I/7. The radiation flux produced at the location will then decreased by exp[—T.y^{rT, e)], where T^-,{rr, e) — 
f°° dzu f°° aryj{-cu, e)n^{r, vj)dr is the photon-photon optical depth and a-,^i{vj, e) is the corresponding cross-section. As ^^^(ro, e) is 

peaked at the threshold energy (Svensson 1987), the above expression can be simplified to T-y.y(rr,e) — {aT/5) f e^i^n^{r,e^.y)dr, 
where n^{r, e^.^) is the photon density at threshold at the location r (see equation|6). 

As ^ 7, this absorption mechanism would be important as long as the companion star produces a sufficient nuinber of photons with 
energies e-^^, > I/7. This limit is illustrated in Fig.|2|for various evolutionary histories of the stellar companion, which has been assumed 
to be in the main sequence (Izzard et al. 2003). It can be seen that the number of soft photons able to interact with the high-energy 7-rays to 
produce pairs strongly increases as the bulk Lorentz factor of the relativistic wind exceeds 10^. The absorbed radiation will subsequently 
be reprocessed by the pairs and redistributed in energy. Each electron and positron will have an energy 7± ~ e/2 at birth, and will cool 
as a consequence of the Compton-drag process. The positrons will in turn annihilate in collisions with electrons in the wind, producing a 
blueshifted annihilation line at e ~ 7. 

Direct measurements of the characteristics of this hard-energy radiation are frustrated by the fact that in traversing intergalactic distances, 
7-rays may be absorbed by photon-photon pair production on the background field radiation (Gould & Schreder 1967). Photons of energy 
near 1 TeV interacting with background photons of 0.5 eV have the highest cross section, although a broad range of optical-infrared 
wavelengths can be important absorbers because the cross section for pair production is rather broad in energy and, in addition, spectral 
features in the extragalactic background density can make certain wave bands more important than the cross section alone would indicate 
(Biller et al. 1998). The current generation of ground-based, 7-ray telescopes have a typical lower energy threshold of ~ 0.5 TeV and are 
thus expected to be able to see sources possessing redshifts up to z = 0.1. The next generation of instruments (e.g. Glast^) is expected to 
have an energy threshold in the region 0.05-0.1 TeV, and will therefore be able to see out to a redshift of at least z = 0.5. Fig.|2|and its 
caption summarise the above limits along with the spectral attributes of the observed 7-ray signal. 
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Figure 4. Luminosity of the scattered 7-ray radiation produced by tlie interaction of the postburst relativistic outflow with the supernova photon field. The 
properties of the observed signal are linked to the postexplosion expansion history of the SN envelope. The bolometric luminosity of SN1998bw, as derived 
by Woosley et al. (1999), is used to calculated the photon field energy and column density (inset panel). We assume v = 10^ cm s~^ and i/> = 0.7. The thick 
dashed line denotes the threshold energy at which the beamed photons are absorbed by the radiation from the SNR. 

3 A DECAYING MAGNETAR IN AN ASYMMETRIC SNR 

The success of the previous model is inseparably linked with the assumed residence of a companion star. Here we consider a related and less 
restrictive scenario in which the ambient photons are produced by the postexplosion expansion of the disrupted envelope. Photons emitted 
from the supernova remnant (SNR) shell provide copious targets for the relativistic outflow to interact and produce high energy 7-rays. 

GRBs are thought to be produced when the evolved core of a massive star collapses either to a fast-spinning neutron star or a newly 
formed black hole. In the latter case, a GRB is likely to be triggered if the remaining star has sufficient angular momentum to form a 
centrifugally supported disk (e.g. MacFadyen & Woosley 1999). A funnel along the rotation axis would have been blasted open during the 
1-100 s duration of the original burst; it would subsequently enlarge owing to the postexplosion expansion of the envelope of the progenitor 
star (e.g. Woosley 1993). The ram pressure of the continuing MHD outflow would further enlarge the funnel. Besides asymmetrically ejecting 
the SNR envelope (e.g., no ejecta in the polar direction), the supernova explosion may leave behind parts of the He core which take longer 
in falling back. Additional target photons may arise from parts of the disrupted He core, no longer in hydrostatic equilibrium and moving 
outwards inside the SNR shell. For a nominal subrelativistic shell speed v = l(fvg cm s~^, the typical distance reached is rsnr ~ W'^^vgtd 
cm in days. The outflowing wind from the compact remnant (which we assume to be relativistic) would propagate inside the funnel cavity 
(of conical shape with semi-aperture angle ip) while interacting with the SNR target photons with typical energy Ssnr ~ 10^^ — 10^^ (i.e. 
optical frequencies). 

Under the foregoing conditions the particles in the postburst relativistic outflow see blue-shifted photons pouring in from the forward 
direction. The rate of energy loss of a relativistic particle moving in a radiation with an energy density Wsnr ~ Lsut / r^^-ivC is about 
nieC? d'y / dt ~ —WsniaTcy^ in the Thompson limit. The total luminosity of scattered hard photons Lr is equal to the total particle energy 
losses in the course of motion from the magnetar to infinity: LriOsm) ~ [A'y{6snT) /j]Ln- Lr varies with 6sni which in turn evolves with 
the expansion history of the stellar envelope. Fig.|4|shows the total luminosity of scattered hard photons Lr as a function of the age of the 
remnant. The bolometric luminosity of SN 1998bw, as derived by Woosley et al. (1999), is used here to calculated the photon field energy 
by assuming a constant expansion velocity of v — 10^ cm s~^. The resulting radiation pressure on electrons in the ejecta will brake any 
outflow whose initial Lorentz factor exceeds some critical value 7cd < {La/Lsnv)^^^, converting the excess kinetic energy into a directed 
beamed of scattered photons. In reality, the external parts of the postburst relativistic outflow, which are in closer contact with the funnel 
walls, are dragged more efficiently since the soft photons arising from the walls can penetrate only a small fraction of the funnel before 
being upscattered by relativistic electrons. The outflow itself, is then likely to develop a velocity profile with higher Lorentz factor along the 
symmetry axis, gradually decreasing as the polar angle increases. Moreover, the outflow power may fluctuate; so also may its baryon content, 
due to entrainment, or to unsteadiness in the acceleration process at the base of the outflow. In this case, additional 7-rays can be produced 
in relativistic shocks that developed when fast material overtakes slower material (i.e. internal shocks). 



4 DISCUSSION 

The initial, energetic portion of the relativistic jet, with a typical burst duration of 10 s, will rapidly expand beyond the envelope of the 
progenitor star, leading in the usual fashion to shocks and a decelerating blast wave. A continually decreasing fraction of energy may 
continue being emitted for periods of days or longer, and its reprocessing by the soft photon field radiation can yield a continuum luminosity 
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extending into the 7-ray band. Photons from a binary companion and/or photons entrapped from a supernova explosion provide ample targets 
for the wind to interact and produce high energy 7-rays. With a sensitivity of 5 x 10~^ ph cm~^ s~^ (6 x 10^^ ph cm~^ s~^) in a '-^ 10^ s 
exposure, the 7-ray instrument on board of Integral (Glast) will detect a 10*^ erg s~^ signal out to a distance of z ~ 0.5 (~ 1). A magnetar 
with Po = 0.8 ms and B ~ 10^'' G would lead to Lsi ~ lO*'^ erg s~^ after 1 day. This could also be a consequence of a drop in B from 
10^'' to 3 X 10^^ G in a compact structure with stored energy of at least 10^^ erg whose characteristic spin period remained constant (at a 
fraction of a millisecond). A similar argument can be developed based on the concept of a viscosity. For a hot dense torus around a black hole 
resulting from collapse of the core of the progenitor star, the viscous accretion time for a torus of radius 10^ rg cm is tvis ~ 1.5rg^^ B^2 5 
for nT ~ constant, and the accretion of < 1O^^M0 in t ~ 1 day is sufficient to provide a characteristic L ~ 10'*^ erg s~^. This luminosity 
may not dominate the continuum GeV afterglow if it is only modestly reprocessed (i.e Lp < 10** erg s^^). In this case, we cannot rule out 
the possibility that much of the GeV emission is due to inverse Compton (Meszaros & Rees 1994) or synchrotron self-Compton (Derishev et 
al. 2001) losses in the standard decelerating blast wave, which could produce an afterglow luminosity that was still, 1 day after the original 
explosion, as high as 10*^ erg s~^ (see e.g. Dermer et al. 2000 for a blast wave expanding into a medium with no = 100 cm~^). 

In closing, the scattered 7-rays from a peculiar Ib/c SNR of greater than usual brightness or a high luminosity stellar companion (i.e. 
> IOM0) should be easily detectable (at tobs < 1-5 days) by the recently launched Integral. For a less luminous accompanying star (i.e. 
< 5Mq), the expected 7-ray flux can be detected up to 2; = 0.1, although may be difficult to disentangle from the scattered SNR emission 
whose intensity is likely to dominate at t < 100 days. The 7-ray signals for "mean" events (i.e. z ~ 1) would stand out with high statistical 
significance above the background provided that tobs < 0.1 — 0.5 days. The planned Glast experiment may also soon provide relevant limits 
(i.e. z < 0.5) for individual events at higher threshold energies (Figs. 2 and 3). Its detection would unveiled the presence of a postburst 
relativistic wind, and the spectral signatures of the 7-ray emission would help constrain the nature of the precursor star. 
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